<n: 

>v 

a ■ 



^o 
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We show that quantum-interference phenomena can be realized for the dissipative nonlinear sys- 
tems exhibiting hysteresis-cycle behavior and quantum chaos. Such results are obtained for a driven 
dissipative nonlinear oscillator with time-dependent parameters and take place for the regimes of 
long time intervals exceeding dissipation time and for macroscopic levels of oscillatory excitation 
numbers. Two schemas of time modulation: (i) periodic variation of the strength of the x(3) nonlin- 
earity; (ii) periodic modulation of the amplitude of the driving force, are considered. These effects 
are obtained within the framework of phase-space quantum distributions. It is demonstrated that 
the Wigner functions of oscillatory mode in both bistable and chaotic regimes acquire negative 
values and interference patterns in parts of phase-space due to appropriately time-modulation of 
the oscillatory nonlinear dynamics. It is also shown that the time-modulation of the oscillatory 
parameters essentially improves the degree of sub-Poissonian statistics of excitation numbers. 

PACS numbers: 42.50.Dv 
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I. INTRODUCTION 

In recent years the study of quantum dynamics of os- 
cillators with time-dependent parameters has been focus 
of considerable attention. This interest is justified by 
many applications in different contexts. Particulary, one 
application concerns to the center of mass motion of a 
laser cooled and trapped ion in a Paul trap [l[. The 
quantum dynamics of an anharmonic oscillator (AHO) 
with time dependent modulation of its frequency and 
nonlincarity parameters has been investigated in appli- 
cations to macroscopic superposition of quantum states 
0. In the last few years there has been rapid progress 
in the construction and manipulation of nanomechanical 
oscillators with giant x(3)-Kerr nonlinearity Q-0. The 
nanomechanical resonator with a significant fourth-order 
nonlinearity in the elastic potential energy has been ex- 
perimentally demonstrated jg|. It has also been shown 
that this system is dynamically equivalent to the Duff- 
ing oscillator with varied driving force [7fl. This scheme 
is widely employed for a large variety of applications as 
well as the other schemes of micro- and nanomechani- 
cal oscillators, more commonly as sensors or actuators in 
integrated electrical, optical, and optoelectrical systems 

It is well assessed that in the case of unitary dynam- 
ics, without any losses, an anharmonic oscillator leads to 
sub-Poissonian statistics of oscillatory excitation number, 
quadratic squeezing and superposition of macroscopically 
distinguishable coherent states. For dissipative dynam- 
ics the important parameter responsible for production 
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of nonclassical states via x(3) materials is the ratio be- 
tween nonlinearity and damping. Therefore, the prac- 
tical realization of such quantum effects requires a high 
nonlinearity with respect to dissipation. In this direction 
the largest nonlinear interaction was proposed in many 
papers, particularly, in terms of electromagnetically in- 
duced transparency [9[ and by using the Purcell effect 
[lOj . and in cavity QED [11]. The significant nonlinear- 
ity has also been observed for nanomechanical resonators 
[12| . These methods can lead to x(3) nonlinearity of 
several orders of magnitude higher than natural optical 
self-Kerr interactions. Note, that high x(3) nonlinear os- 
cillators generate also a lot of interest recently due to 
their applications in areas of quantum computing [13]. 

In the case of nonlinear dissipative x(3) interaction 
stimulated by coherent driving force, the time evolu- 
tion cannot be solved analytically for arbitrary evolu- 
tion times and suitable numerical methods have to be 
used. Nevertheless, with dissipation included a driven 
AHO model has been solved exactly in the steady-state 
regime in terms of the Fokker-Planck equation in complex 
P representation [14j. Analogous solution has been ob- 
tained for a combined driven parametric oscillator with 
Kerr nonlinearity [l5( . The Wigner functions for both 
these models have been obtained using these solutions 

MM- 

The investigation of quantum dynamics of a driven dis- 
sipative nonlinear oscillator for non-stationary cases is 
much more complicated and only a few papers have been 
done in this field up to now. More recently, the quan- 
tum version of dissipative AHO or the Duffing oscillator 
with time-modulated drivin g fo rce has been studied in 
the series of the papers [l7| , [l8| , [l9| in the context of a 
stochastic resonance [l7| , quantum-to-classical transition 
and investigation of quantum dissipative chaos |18| , [191 ] . 

In this paper we continue investigation of time- 
modulated effects for quantum version of the Duffing os- 



cillator. In this direction we investigate quantum effects 
in the presence of dissipation and decoherence, mainly 
sub-Poissonian statistics and signature of quantum inter- 
ference within the framework of the Wigner function for 
an oscillatory mode. Two schemas of time modulation: 
(i) periodic variation of the strength of the x(3) nonlin- 
earity; (ii) modulation of the amplitude of the driving 
force will be considered. Our main result is that time 
modulation of the oscillatory parameters essentially im- 
proves the degree of sub-Poissonian statistics of oscilla- 
tory excitation numbers as well as leads to negative val- 
ues of the Wigner functions in an over transient regime 
for definite time intervals exceeding the transient time 
within the period of modulation. Thus, we demonstrate 
that for the case of time-modulated AHO the Wigner 
function gradually deviates from the Wigner function 
corresponding to an ordinary AHO without any mod- 
ulation. Really, surprisingly simple analytical results for 
the Wigner functions of an ordinary AHO have been ana- 
lytically obtained in over transient regime [14| , [16| , that 
are positive in all ranges of phase-space. In our scheme 
the negative values for the Wigner function of oscilla- 
tory mode, that reflect quantum-interference patterns in 
phase space, appear due to time-modulation dynamics 
and are realized in both bistable and chaotic operational 
regimes. 



Note, that nowadays, large theoretical and experimen- 
tal papers have been devoted to the implementation of 
macroscopic (i.e., many particle) quantum superposition 
states (MQS). The most notable results on MQS are ex- 
perimentally obtained with atoms interacting with mi- 
crowave fields trapped inside a cavity [2(| or for freely 
propagating fields [21( and with optical parametric am- 
plifier in the presence of decoherence [22[ • Naturally, the 
question arises whether the quantum-interference pat- 
terns obtained for driven AHO can be realized for time in- 
tervals exceeding the characteristic decoherence and dis- 
sipation time on many particle level. The correspond- 
ing results obtained here are verified for the operational 
regimes of high nonlincarity with respect to dissipation, 
e.g., for the ratio xll — 1 ' 10 _1 "=" 1-5 • 10~ 3 , where \ 
is the nonlinear constant proportional to x(3) and 7 is 
the damping constant. Nevertheless, we have checked for 
these parameters that phase-space interference patterns 
are also obtained in the macroscopic level that involves a 
number of oscillatory excitation numbers in excess of 52. 



The outline of this paper is as follows. In the next 
section we describe the models and investigate the mean 
oscillatory excitation numbers and quantum fluctuations 
of excitation numbers on the base of the Mandel param- 
eters. In Sec. Ill we shortly discuss the case of the uni- 
tary dynamics. In Sec. IV we present the results for 
the Wigner functions showing negativity due to the time 
modulations of the oscillatory dynamics. We summarize 
our results in Sec. V. 



II. MODELS: EXCITATION NUMBERS AND 
QUANTUM STATISTICS 

We treat the Duffing oscillator as an open quantum 
system and assume that its time evolution is described 
by Markovian dynamics in terms of the Lindblad master 
equation for the reduced density matrix p. In the in- 
teraction picture that corresponds to the transformation 
p — > e -iua at p e iuia at ^ wnere a + anc [ a are t ne Bose an- 
nihilation and creation operators of the oscillatory mode 
and tu is the driving frequency, this equation reads in the 
Markovian form as 

dp —i , 



E ( L *P L t - \ L tLiP - \pLiL^) . (1) 



i=l,2 



The Hamiltonians are 



H = HAa + a, 

H mt = h X (t)(a+a) 2 + h(f(ty 
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(2) 



where x(t) an d f(t), which may or may not depend on 
time, represent, respectively, the strength of the nonlin- 
earity and amplitude of the force, uj is the resonant fre- 
quency, A = luq — uj is the detuning. The dissipative 
and decoherence effects, losses, and thermal noise are in- 
cluded in the last part of this equation, where Li are the 
Lindblad operators: 



Li = y/(N + l)ia, 



L 2 = 



(3) 



7 is the spontaneous decay rate of the dissipation process 
and N denotes the mean number of quanta of a heat bath. 

This model seems experimentally feasible and can be 
realized in several experimental schemes. In fact, a sin- 
gle mode field is well described in terms of an AHO, and 
the nonlinear medium could be an optical fiber or a x(3) 
crystal, placed in a cavity. The anharmonicity of mode 
dynamics comes from the self-phase modulation due to 
the photon-photon interaction in the x(3) medium. In 
this case, it is possible to realize time modulation of the 
strength of the nonlinearity by using a media with peri- 
odic variation of the x(3) susceptibility. 

On the other side, the Hamiltonian described by Eq. 
@ describes a single nanomechanical resonator with a + 
and a raising and lowering operators related to the posi- 
tion and momentum operators of a mode quantum mo- 
tion 



2mu>o 



(a + a + ), p = — zy 2hmuio(a — a ), (4) 



where m is the effective mass of the nanomechanical res- 
onator, ujq is the linear resonator frequency and x propor- 
tional to the Duffing nonlinearity. One of the variants of 
nano-oscillators is based on a double-clamped platinum 
beam [12J for which the nonlinearity parameter equals to 



X = h/4:\/?>Qma 2 cl where a c is the critical amplitude at 
which the resonance amplitude has an infinite slope as 
a function of the driving frequency, Q is the mechani- 
cal quality factor of the resonator. In this case, the gi- 
ant nonlinearity \ — 3.4 • 10 _4 s _1 was realized. Note, 
that details of this resonator, including expression for 
the parameter a c , are presented in |23|. On decreasing 
nanomechanical resonator mass, its resonance frequency 
increases, exceeding 1 GHz in recent experiments Q, [Ji- 
lt is possible to reach a quantum regime for such fre- 
quencies, i.e., to cool down the temperatures for which 
thermal energy will be comparable to the energy of oscil- 
latory quanta. The recent investigations in this direction 
are devoted to classical to quantum transition of a driven 
nanomechanical oscillator [24j , generation of Fock states 
25], nonlinear dynamics, and stochastic resonance |26j . 

Cyclotron oscillations of a single electron in a Penning 
trap with a magnetic field are another realization of the 
quantum version of the Duffing oscillator [27| , [28| , [29j ■ 
In this case the anharmonicity comes from nonlinear ef- 
fect that is caused by the relativistic motion of an elec- 
tron in a trap, while the dissipation effects arise from the 
spontaneous emission of the synchrotron radiation and 
thermal fluctuations of the cyclotron motion. Note that 
a one-electron oscillator allows one to achieve a relatively 
strong cubic nonlinearity, x/l ^5 1- 

For the constant parameters %(£) = \ an d fit) = f 
the equations ([T]) and ([2j describe the model of a driven 
dissipative AHO that was introduced long ago in quan- 
tum optics to describe bistability due to a Kerr nonlinear 
medium [30J. For the case of time-dependent parame- 
ters x(t) and f(i) the dynamics of the AHO exhibits a 
rich phase-space structure, including regimes of regular, 
bistable and chaotic motion. We perform our calculations 
for these three qualitatively different regimes concern- 
ing two models of time-modulated AHO corresponding 
to two physical situations: (i) \ = x(t) — Xo + x.i s ^ n (^t) 
and f(t) = const = /; (ii) / = f(t) = f + /ism(ftt) and 
x{t) = const = x with S -C uj and $7 <C u the modulation 
frequencies. 

We analyze the master equation numerically using 
quantum state diffusion method (QSD) [3l|]. According 
to this method, the reduced density operator is calculated 
as the ensemble mean 
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P (t) = M(|^(t))(^(t)|) = lira - ^ IMQWM 

(5) 
over the stochastic pure states \4>t(t)) describing evolu- 
tion along a quantum trajectory. The stochastic equation 
for the state \^(t)) involves both Hamiltonian described 
by Eq. ([2]) and the Linblad operators described by Eq. 
Q. We calculate the density operator using an expan- 
sion of the state vector \i/j^) in a truncated basis of Fock's 
number states of a harmonic oscillator 



hfc(*)) = £a£(*)|n). 




FIG. 1: The photon excitation numbers for two cases: (1) 
for monostable regime with parameters A/7 = 0.1,x(t)/7 = 
5 • 1CT 3 (1 + 0.7sin(<5t)),d/7 = 3,f/j = 20; (2) for bistable 



regime with parameters A/7 = 
0.5sin(<5i)),<5/7 = 3,// 7 = 20. 



-i,x(t)h = 6-10- 3 (l + 



A. Improvement of the sub-Poissonian statistics by 
modulation of the nonlinearity 

In the classical limit the corresponding equation of mo- 
tion for the dimensionless mean amplitude a = Tr(pa(t)) 
has the form 



da 

It 



7 



■i(A+( Xo + Xisin^)(l + 2|a| 2 ))a- J /. (7) 



(6) 



We remind the reader that for the constant strength of 
the nonlinearity, xi = 0\ the semiclassical mean excita- 
tion number \a\ 2 as a function of the driving intensity 
can exhibit multiple steady states and hysteresis. In the 
case of modulation, xi 7^ 0, the interaction Hamiltonian 
is explicitly time dependent, and the system exhibits re- 
gions of regular, bistable, and chaotic motion with xoj 
Xi, A, and the modulation frequency being the control 
parameters. 

We examine these operational regimes by numerical 
analysis of the phase space of dimensionless position and 
momentum x — Re(a) and y = Im(a). Choosing xq 
and yo as an arbitrary initial phase-space point of the 
system at the time to, we define a constant phase map 
in the (X, Y) plane by the sequence of points (X n ,Y n ) = 
(X(t n ),Y{t n )) at t n = t + ^-n (n = 0,1,2,...). This 
means that for any t — t n the system is at one of the 
points of the Poincare section. The analyses show that for 
time scales exceeding the damping time I/7, the asymp- 
totic dynamics of the system is regular in the limits of 
small and large values of the modulation frequency, i.e., 
i « 7, 1$ > 7, and for positive values of the detuning, 
A > 0. Such regime is also realized when the modulation 
part of the parameter, i.e., Xi is much less than the xo- 

First, we discuss the case of classically regular be- 




Comparing this result with the case without any time- 
modulation we also present the Q parameter for the max- 
imal and minimal values of x(t)i he., for x — 8-5 ■ 10~ 3 
(Q = -0.3, curve (b)) and for x = 1.5- 1CT 3 (Q = -0.36, 
curve (c)). Thus, we conclude that the modulation of 
the strength of nonlinearity leads to the improvement of 
the level of sub-Poissonian statistics below the level of 
corresponding steady-state regime. Indeed, we find that 
Qmin < —0.36. It is interesting that Qmin corresponds 
to the maximal value of the mean excitation number, 
(n) = 249. As calculations show, analogous result takes 
place for the bistable regime, Fig. [TJ case(2). In this case, 
Qmin — —0.7 for the mean excitation numbers (n) = 352. 
Note that such effect of improving the sub-Poissonian 
statistics has been recently obtained for the dissipative 
AHO under time-modulation of the driving amplitude 
181. 



FIG. 2: The Mandel parameter for the monostable regime 
with the parameters as in Fig. [JJ case (1), (curve(a)); AHO 
without time-modulation: \H = 8-5 ■ 10 -3 (b); \ll = 1-5 ■ 
10 (c). The other parameters are as in the Fig. [TJ case (1). 



III. TIME MODULATION IN THE CASE OF 
UNITARY DYNAMIC 



havior assuming the interaction of the system with vac- 
uum reservoir: N = 0. Below we present numerical 
results for both the oscillatory mean excitation number 
(n) = A! ((ip^\a + a\tp^)) and Mandel Q parameter which 
describes the deviation of excitation number uncertainty 
from the Poissonian variance, i.e., Q = — — "L\ , 



((An) 2 ) - ((a+a) 2 ) - (a+a) 2 , i.e., 

((An) 2 )=M({^\(a + a) 2 m)-{n) 



(8) 



For the driven dissipative AHO under time-modulation 
of the strength of nonlinearity, the ensemble-averaged 
mean oscillatory excitation number exhibits a periodic 
time-dependent behavior in both cases of regular and 
chaotic regimes for over transient time-intervals. For 
bistable and monostable regimes the results of numer- 
ical calculations for time evolution of (n) are depicted 
in Fig. [TJcase 1) and Fig. [TJcase 2) for the parame- 
ters: %{t)h = 5 ' 10~ 3 (1 + 0.7sin(37t)) and x(t)/ 7 = 
6 ■ 10~ 3 (1 + 0.5 sin(37t)). The time evolution of the Man- 
del parameter is depicted on Fig. [2j As we see (Fig. 
[2Ka)), the Q parameter also shows a time-dependent 
modulation and formation of the sub-Poissonian statis- 
tics (Q < 0) for the definite time intervals exceeding 
the transient time by > 6. The level of the oscilla- 
tory excitation-number fluctuations reaches to the min- 
imum values Qmin = —0.69 and the maximum values 
Qmax = 1-2 at the fixed intervals of time, respectively, at 
jt k = 6.8 + 2^7 and at jt k = 7.6 + ^7 , {k = 0, 1,2...). 



We now calculate the quantum distributions, at first 
considering the simplest case of non-dissipative AHO 
without any driving. In this case the operator (a + a) 
is an invariant of motion, since this operator commutes 
with the corresponding Hamiltonian H(t) — h[Aa + a + 
x(t)(a + a) 2 ]. The system exhibits an unitary evolution; 
for an initial coherent state \ao) with the phase <p = at 
t = 0, the state evolution at time t is 



I^(t)) = j7(t)|oo; 



- p-M/2 



E 

ra=0 



"0 



^W" 2 |n)» (9) 



, where tp(t) = J x(j)dr. This state involves the time- 
modulation effects as well as describes the interference of 
the Fock number states \n). The corresponding nonsta- 
tionary Wigner function is written as 

We transform this formula to the form 

W(r,6) = Y,Pnm(t)W mn (r,8), (11) 



where (r, 9) are the polar coordinates in the complex 
phase-space plane, x = rcosO, y — rsin9, while the co- 
efficients W mn (r, 9) are the Fourier transform of matrix 
elements of the Wigner characteristic function 



W mn (r,0) = 




(4r ), n > m 



(12) 



r 



The matrix elements p n 
equal to 



,(t) - (r#(t))<V(*)W 



I I ti-\-tti 
O (t) - e -| Q "| 2 ™ e iv(t)(n 2 -m 2 ) 

V n!m! 
where for x(i) = xo + Xisin(St + (f>) we have 



¥>(*) =X0< + -y[c08((5t 



COS( 



(13) 



(14) 



These results arc valid for short time intervals, much 
less than the characteristic dissipation time, t -C 7 -1 . 
The Wigner function Eq. (jlljl with the density matrix 
Eq. (|13|) leads to the interference fringes arising from 
the non-diagonal elements p nm - For non- modulated case, 
f{t) = Xo^the Wigner function is well studied (see, for 
example, [32]). Particularly, for ~xt=\ and in the case of 
the initial coherent state, it describes a superposition of 
other coherent states of the same amplitude but different 
phases. The novelty of the Eqs. (|11|) and (fT3|) consists in 
the time-dependent nonlinearity that is the factor of f{t). 
This modulation term allows to control the time intervals 
of the maximal superposition with appropriate choice of 
the modulation frequency and the phase </>. Indeed, the 
first time at which we have the superposition of coher- 
ent states is determined from ip(t) = £. Particularly, for 
St -C 1, we have ip(t) ~ (xo + Xi sm(j))t and hence the 
superposition is realized, if t = ir/2(xo + Xi sin</>). 



IV. WIGNER FUNCTIONS AND 
QUANTUM-INTERFERENCE PATTERNS 

It is well known that the phase-space Wigner distri- 
bution function can simply visualize nonclassical effects 
including quantum-interference. For example, a signa- 
ture of quantum interference is exhibited in the Wigner 
function by the non-positive values. In this section the 
numerical results of the nonstationary Wigner functions 
in bistable and chaotic regimes of AHO are presented and 
discussed. 

It should be noted that the most of investigations of 
the quantum distributions of oscillatory states, including 
also modes of radiations, have been made for the steady- 
state situations. The simplicity of Kerr nonlinearity al- 
lows to determine the Wigner function of the quantum 
state under time evolution due to interaction (see, for ex- 
ample [32]). In this sense, we note the main peculiarity 
of our paper in comparison with above noted important 
inputs. In this paper, we calculate the Wigner functions 
in an over transient regime, t 3> 7 -1 , of the dissipative 



dynamics, however, we consider time-dependent effects 
which appear due to the time-modulation of the oscilla- 
tory parameters. 

Below we investigate the Wigner functions for both 
cases of time-modulation (see, cases (i) and (ii)). For 
periodic variation of the strength of x(3) nonlinearity, 
scmiclassical dynamics is described by Eq. ([7]), while for 
the case of time-modulation of the driving amplitude the 
scmiclassical equation reads as 



da 

~~dt 



7 
--a 



i(A + X o(l + 2|a| 2 ))a- 



i(fo + fisin(Qt)). 

(15) 



A. Quantum interference pattern assisted by the 
bistability. 

First, we consider the case of time-modulated driving 
force. In the limit /1 <C /o the system is reduced effec- 
tively to the model of a single driven anharmonic oscil- 
lator, which exhibits bistability for the definite range of 
the parameters: xo, A, /o, and 7 (see for example [7|, 
[17J |. [29j j. [301), however, for negative detuning, A < 0. 
In this case, the semiclassical steady-state oscillatory ex- 
citation number |a| 2 , as solution of Eq. (|j"5l) . displays 
hysteresis behavior, while quantum mechanical mean ex- 
citation number (a + a) docs not show any hysteresis in 
the critical range due to quantum statistical averaging 
fl4| . The analysis of the stochastic trajectories for an 
expectation number n^(t) = (-0^|a + a|-0^) shows that the 
system spends most of its time close to one of the semi- 
classical bistable solutions with quantum-interstate tran- 
sitions, occurring at random intervals [17| . On increasing 
the amplitude /1 a new channel of stimulated interstate 
transitions is raised between the semiclassical bistable so- 
lutions. Their contribution for the further increasing of 
/1 at /1 < /o leads to the emergence of a chaotic regime. 

In another scenario of transition from regular regime 
to bistability and chaos the modulation frequency fl may 
be varied, with the other parameters unchanged. In the 
range Q ^C 7 the modulation of the system is adiabatic. 
On increasing modulation frequency the system oscillates 
between the two possible metastable states. At f2 > 7 a 
strong entanglement of these states occurs, and the sys- 
tem comes to a chaotic regime. As can be seen from the 
numerical calculations, the analogous results are also ob- 
tained from the case (i), i.e., for the case of modulation 
of the parameter x(i). Some differences between oper- 
ational regimes of these two models will be considered 
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FIG. 3: The Wigner function and its contour plot, averaged 
over 3000 trajectories, for the parameters: A/7 = —15, xll = 
2, //7 = 5.8(1 + 0.5 sin(fit)), fi/7 = 2, (a) and (b). The con- 
tour plot of the Wigner function that indicates the ranges of 
negativity in the black (b). The Wigner function for AHO 
without modulation, //7 = 5.8 (c). The probability distribu- 
tion of excitation numbers (d): case of time-modulation (1), 
case without of modulation (2). The parameters are as on the 
Figs, (a) and (c), respectively. 



FIG. 4: The Wigner function (a); the contour plot of the 
Wigner function (b); and the mean excitation number (c); 
averaged over 3000 trajectories, for the dissipative AHO un- 
der time-modulated driving. The parameters are in the 
range of bistability: A/7 = -13.02, xll = °- 08 > fll = 
29(1 + 0.5sin(fi£)),a/7 = 2. The hysteresis curve (d) for 
the stationary case. 
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FIG. 5: The Wigner function (a); the contour plot of 
the Wigner function (b); and the distribution of quadra- 
ture amplitude P(x) for the driven dissipative AHO with 
time-modulation of the nonlinearity. The parameters are 
in the range of bistability: A/7 = — 15, x(t)/l — 1-5(1 + 
0.5sin(5i)),<5/7=5,//7 = 10.2. 



below. 

The Wigner function is determined by the expression 



W(a,t) 



A e -2M 2 

7T 2 



d 2 p{-(i\p{t)\fi)e 



-2(j3a*-l3*a) 



( 16 ) 

which is obtained from the Eq. dill) , where, however, the 
matrix elements of the density operator of the full dissi- 
pative system are used, pnm{t) = M((n\ip^(t))(ip^(t)\m}). 



The coefficients W mn (r, 6) are the Fourier transform of 
matrix elements of the Wigner characteristic function 
(see Eq. (fT2jl). 

As calculations show, the most striking signature of 
quantum bistability in the presence of time-modulation 
is the appearance of quantum-interference with negative 
regions in the Wigner function. We illustrate quantum- 
interference pattern in phase space and in the bistable 
operational regime on Figs. [SJIU More importantly, we 
clearly observe (Fig. Efa)) that the Wigner function dis- 
plays two peaks around the bistable semiclassical solu- 
tions as well as the negative part between them that in- 
dicates quantum-interference. The ranges of negativity 
are indicated in details on Fig. EJb). Another pecu- 
liarity is that the Wigner function is nonstationary and 
the interference patterns take place for the definite time 
intervals jtk = 6.9 + ^f^7 (k=0,l,2...), exceeding tran- 
sient time, at which the mean excitation number (n) 
reaches its maximal value. The interference pattern is 
destroyed as the time-modulation is decreased. Indeed, 
as it is shown on Fig. (He), in the case /1 = the Wigner 
function is positive in all phase-space and indicates quan- 
tum bistability in accordance with the analytical results 
(14| . Examples of the curves for the ensemble-averaged 
excitation number distributions P n = (n\p\n) are demon- 
strated on Fig. [3jd) for the two cases of modulated and 
non-modulated dynamics described on the Figs. 3(a) 
and 3(c). In correspondence with the bistable behav- 
ior of the system we find a developing bimodal struc- 
ture of P n — (n\p\n) for the oscillatory mode. However, 
this structure is less pronounced for the case of non- 
modulated stationary dynamics and is well resolved for 
the case of time-modulation of the driving force, as we 
can also observe from the results for the Wigner func- 
tions. The reason is that the interstate transitions due 
to time-modulation of the driving force approximately 
equalizes the populations of the states. 

Note, that these results are given for a choice of pa- 
rameters which ensures that the oscillator is in the spe- 
cific bistable regime operating with small main excita- 
tion numbers. Nevertheless, it is possible to find the 
parameters that would enable us to observe quantum- 
interference pattern on a macroscopic level of excitation 
numbers. For this goal, we use the scaling properties of 
the driven AHO recently proposed in [18|]. Indeed, it is 
easy to verify that Eqs. (O and (|15l) are invariant with 
respect to the following scaling transformation of the am- 
plitude a — > a — Aa, where A is a real positive dimen- 
sionless coefficient, if the oscillatory parameters are cor- 
respondingly transformed as A ^ A = A + x(l — 1/A 2 ), 

X -> X = XA 2 , / -> / = A/, 7 -> 7' = 7- This 
scaling properties of the semiclassical equations signifies 
that for the definite different sets of the parameters the 
phase-space trajectories have the same form and differ 
from each other only in a scale. It has been numerically 
shown that this scaling parameter is also approximately 
realized in a quantum ensemble theory in the presence 
of quantum noise, for wider ranges of the parameters, 



but not for large values of the parameter x/7, i.e., in a 
deeper quantum regime. We calculate the Wigner func- 
tion and its contour plot for the parameters obtained 
from the data of the Fig. 3(a) with the scaling parameter 
A = 5. The results presented in Fig. BJa) and Fig. HJb), 
particularly, illustrate the quantum interference pattern 
for a macroscopic level of the mean excitation number, 
(n) = 40. As we see, the negative part of the Wigner 
function decreases with increasing the mean oscillatory 
excitation number. The result for time evolution of the 
mean oscillatory excitation number averaged over quan- 
tum trajectories is depicted in Fig. BJc). For such pa- 
rameters the system still operates in a bistable regime. 
We indicate the hysteresis curve for the corresponding 
semiclassical solution for the excitation number \a\ 2 on 
Fig. Bid) for the stationary case: /1 = 0. 

Analogous results are obtained for the case (i), that 
is the time-modulation of the strength of the nonlinear- 
ity. The results of numerical calculations of the Wigner 
function and the corresponding distribution of quadra- 
ture amplitude P(x) in the operational regime of bista- 
bility are shown on Fig. 01 

Note that these interference effects have been demon- 
strated for the strong quantum regimes of high nonlinear- 
ity. As a consequence, the oscillatory excitation number 
is small in this bistable regime. Below we consider the 
regime of chaotic dynamics where quantum-interference 
pattern is realized with relatively high level of excitation 
numbers. 




Re(a) 



(c) 



B. Negative Wigner functions for the chaotic 
dynamics. 

For both models of AHO corresponding to two forms 
of time- modulations (cases (i) and (ii)), the dynamics 
of the system is chaotic in the ranges: S > 7 or Q > 
7 and xo — Xi ° r /o — /1 an d for negative detuning. 
Thus, the ways to realize the controlling transition from 
the regular to chaotic dynamics through the intermediate 
ranges of bistability are to vary the strength \i or /1 of 
the modulation processes in the ranges from xi ^ Xo 
or /1 « /o to xi < Xo or /1 < / . In order to get 
a quantitative measure for chaotic dynamics, we have 
chosen the analysis on the base of the Poincare section. 

We note, that the order-to-chaos transition for dissi- 
pative AHO under time-modulated driving force was al- 
ready analyzed 1$]. Therefore, here we concentrate on 
investigation of quantum interference phenomenon con- 
sidering AHO with time-modulated nonlinearity. The 
results of the ensemble-averaged numerical calculation of 
the mean excitation number, the Poincare section and 
the Wigner function are shown on Figs. [Ua)-[nic), re- 
spectively. The mean excitation number of the driven 
AHO versus dimensionless time is depicted in Fig. [HJa). 
In contrast to the semiclassical result its quantum en- 
semble counterpart (Fig. Ufa)) has clear regular periodic 
behavior for time intervals exceeding the characteristic 




FIG. 6: (a) The mean excitation number; (b) the Poincare 
section (~ 20000 points) for the dimensionless complex am- 
plitude a, plotted at times of the constant phase Stt = 2nk 
(k=0,l,2...), when the maximal interference pattern on the 
Wigner function (c) is realized, for the case of time-modulated 
nonlinearity xM- The parameters are in the range of chaos: 
A/7 = -5, X{t)/l = 0.2(1 + 0.75 sm(dt)),S/j = 3, // 7 = 10. 



dissipation time, due to ensemble averaging. The Fig. 
IHUb) clearly indicates the classical strong attractors with 
fractal structure that are typical for a chaotic Poincare 
section. Thus, the Wigner function Fig. [SJc) reflects the 
chaotic dynamics, its contour plots in the (x,y) plane are 
similar to the Poincare section. However, the Wigner 
functions have regions of negative values for the definite 
time intervals. The example depicted on Fig. Hfc) corre- 
sponds to time intervals jt/c — 6 + ^f^7 (k=0,l,2...), for 
which the mean excitation number reaches a macroscopic 
level, i.e., n = 52. 



V. CONCLUSION 

We have numerically studied the phenomena at the 
overlap of bistability, chaos and quantum-interference 
for the driven and damped AHO with time-dependent 
parameters. We have pointed out, that the time- 
modulation of the oscillatory parameters, that are the 
strength of x(3) nonlinearity or the amplitude of the driv- 
ing force, leads to formation of the quantum-interference 
patterns in phase-space in over transient regimes, for the 
definite time intervals exceeding the transient dissipation 



time t 3> 7 _1 . These effects are displayed as the negative 
values of the Wigner functions in phase-space. Quan- 
tum interference patterns take place for both bistable 
and chaotic operational regimes and come from the os- 
cillation between possible metastable states of semiclas- 
sical dynamics due to time-modulation. We have also 
demonstrated that the time-modulation of the nonlin- 
earity parameter for the regular operational regime of 
AHO, essentially improves the degree of sub-Poissonian 
statistics of oscillatory mean excitation number. In this 
spirit, we emphasize that the idea of improving the degree 
of quantum effects as well as obtaining qualitatively new 
quantum effects due to appropriately time-modulation of 
open quantum systems was recently exploited for forma- 
tion of high degree continuous variable entanglement in 
nondegenerate optical parametric oscillator [331 ] . 
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